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SUMMARY 


The  question  of  column  end  fixity  is  critically  rec ensile- re ; :.  An  effective 
method  for  determination  is  established.  This  method  relies  upon  the  fact 
that  there  exists  a  readily  determined  point  in  the  span  of  tr.e  eo-urn  lor 
which  the  deflection  due  to  a  given  side  force  is  maximum.  The  product  of 
the  coefficient  which  defines  the  flexibility  at  this  point  and  the  Euler 
load  appropriate  to  the  actual  end  conditions  is,  to  all  intents  and  pur¬ 
poses,  constant. 
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INTRODUCTION 


Discrepansies  between  the  conditions  assumed  in  theory  ani  those  achieved 
in  the  laboratory  or  in  practical  application  seriously  influence  the 
agreement  between  theoretical  predictions  and  actual  observations  on  the 
stability  of  structures.  Admittedly,  methods  of  data  interpretation  which 
enable  valid  deductions  with  regard  to  the  behavior  of  perfect  or  ideal 
structures  to  be  made  from  information  obtained  from  tests  on  imperfect 
or  realistic  structures  have  been  developed.  1*2,  3  Nevertheless,  there 
are  still  some  aspects  of  the  question  that  merit  further  review  and  con¬ 
sideration.  Prime  areas,  in  this  regard,  are  to  be  found  in  the  subjects 
of  boundary  restraint  and  in  the  development  of  nondestructive  testing 
methods  that  can  give  information  with  regard  to  actual  behavior  under 
service  conditions. 

The  subject  of  boundary  restraint  has  long  been  recognized  for  the  simple 
column.  Indeed,  since  the  very  earliest  tests, ^  experimentalists  have,  as 
the  literature  amply  illustrates,  5*  °  striven  to  develop  end  fixtures  or 
restraint  devices  that  would  approximate  as  closely  as  possible  the  bound¬ 
ary  conditions  that  the  analyst  can  treat.  That  they  have,  in  many  cases, 
closely  reached  these  ideal  conditions  is  borne  out  by  the  experimental 
load  values  that  they  have  achieved.  However,  it  is  equally  clear  that  in 
no  case  do  these  special  devices  even  begin  to  approximate  the  conditions 
met  in  real  structures.  We  must  admit  that  theoretical  and  experimentally 
achieved  boundary  conditions  are  often  very  much  apart  from  conditions  of 
realistic  restraint. 

The  problem  is  seriously  magnified  when  the  instability  of  a  plate  is 
examined,  for  even  if  the  plate  is  uniform,  homogeneous,  and  isotropic, 
considerable  difficulty  is  experienced.  This  viewpoint  is  readily  con¬ 
firmed  by  a  review  of  the  current  literature  on  experimental  studies  on 
plates.  When  the  subject  of  sandwich  plates  is  considered,  the  difficul¬ 
ties  are  compounded.  The  conditions  of  edge  restraint  are  much  more 
complex  for  structures  of  this  type  than  for  thin  plates,  a  point  clearly 
made  by  Benson  and  Mayers. 7 

With  regard  to  the  subject  of  shells,  the  matter  is  no  less  important; 
however,  it  has  received  attention  only  in  more  recent  times.  The  majority 
of  the  work  in  this  area  is  analytical,  but  there  are  several  important 
contributions.®*^  10,11  12  jn  stressing  the  importance  of  the  subject 
with  regard  to  this  class  of  structure,  the  earlier  recognition  of  the 
problem  by  Love  and  Southwell^*  1^  cannot  be  ignored. 

In  a  previous  study^  the  various  attempts  that  the  experimentalist  has 
made  to  duplicate,  in  the  laboratory,  the  conditions  prescribed  by  thf 
analyst  have  been  reviewed.  This  work  was  not  intended  to  extend  or  refine 
these  studies.  A  different  approach  was  followed.  From  the  outset,  the 
viewpoint  has  been  that  the  designer  creates  that  which  has  never  been; 
the  analyst  investigates  that  which  has  been  devised.  Thus,  the  axiom  on 
which  the  study  was  based  is  that  the  prime  purpose  of  analysts  mist  be  to 
devise  methods  of  predicting  with  reliability  the  behavior  of  systems  which 


1 


are  subject  to  the  constraints  of  reality.  Mathematical  analysis  and  ex¬ 
perimental  study  are  considered  valid  only  insofar  as  they  lead  to  evalua¬ 
tions  in  this  sense;  in  any  other  fashion  they  are  purely  academic  exercises. 
In  the  literature  that  exists  on  the  influence  of  boundary  conditions  on 
both  stability  and  strength  of  structures,  this  point  is  frequently  lost 
sight  of.  All  too  often  the  theoretician  gives  no  indication  of  hew  his 
prescribed  boundary  conditions  might  be  achieved  in  practice  or  the  degree 
of  error  resulting  from  deviations.  On  the  other  hand.  Idle  experimentalist 
usually  tends  to  conr  ?ntrate,  not  on  determining  the  parameters  which  define 
the  achievable  bouir’  les,  but  rather  on  attempting  to  reproduce  the  condi¬ 
tions  which*  the  th>  cician  lias  specified.  The  practicing  engineer,  the 
designer,  requires  methods,  experimental  and  theoretical,  of  obtaining  a, 
relationship  between  easily  determined  parameters  and  actual  restraint. 

In  this  report,  we  establish  a  rationale. 
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DISCUSSION  OF  TEEOKY 

We  begin  this  study  by  asking  what  is  the  influence  of  boundary  restraint 
on  the  deflection  and/or  stability  of  a  simple  member  under  various  load¬ 
ing  conditions? 


We  note,  in  the  case  of  a  column  under  axial  compression,  that  when  this 
member  is  pinned  at  both  ends,  the  critical  load  level  is  given  by  the 
expression 


P 


cr 


(1) 


provided  the  strut  is  reasonably  slender.  If  the 
changed  from  pinned  to  fixed,  then  the  load  level 
It  becomes  « 

p  _  Vrr  ei 

cr 


end  fixity  condition  is 
is  substantially  altered. 

(2) 


Now,  consider  the  same  simple  member,  but  instead  of  applying  a  destabi¬ 
lizing  compressive  force,  apply  a  lateral  force  at  its  midpoint.  There  is 
no  critical  load  level  for  this  force  system  from  a  stability  viewpoint, 
but  there  is,  of  course,  a  critical  stress.  So  long  as  the  load  is  not 
approaching  the  level  that  would  induce  this  stress,  the  deflection  is 
linear  with  load  and,  for  the  two  end  fixity  cases  considered,  is  given  by 


6  = 


(pinned  ends) 


(3) 


6 


192  EI 


(fixed  ends) 


(M 


It  is  clear  from  these  readily  established  relationships  that  the  change 
in  end  fixity  from  pinned  to  fixed  causes  an  increase  in  critical  load 
level  and  a  decrease  in  deflection  for  a  given  force.  It  is  interesting 
to  note  that  the  destabilizing  force  grows  by  a  factor  of  four;  the 
deflection  decreases  to  one-fourth.  Consequently,  for  these  two  config¬ 
urations,  the  identity 


P 


cr 


(5) 


is  established. 

6 

Now  the  quantity  —  is  in  essence  a  flexibility  coefficient,  and  so  we  may 
say  for  these  two  xases,  at  any  rate,  that  the  product  of  the  critical 
load  for  a  centrally  loaded  slender  column  and  the  flexibility  coefficient 
for  a  central  lateral  force  acting  on  said  column  is  dependent  only  upon 
the  length  of  the  column.  Symbolically: 
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ttl 

TS~ 


--  Constant 


(6) 


P  x’  = 


cr 


The  importance  of  this  relationship,  if  it  can  be  extended  to  other  cases, 
is  clear.  It  provides  the  simplest  possible  method  of  determining  critical 
load  from  easily  obtained  quantities. 


It  is  important  to  note  that  the  tvo  examples  given  so  far  have  certain 
features  in  common: 

1.  They  are  for  geometrically  uniform  bodies. 

2.  They  are  for  bodies  vith  symmetry  and  ideal  end  support. 


To  extend  the  concept  further,  ve  begin  by  systematically  removing  these 
conditions. 

In  the  first  case,  the  restriction  to  nonuniform  bodies  is  removed;  in 
subsequent  sections,  the  support  conditions  are  varied  to  reflect  realistic 
restraint,  including  intermediate  support. 


NONUNIFORM  COLUMNS 


To  begin,  ve  consider  the  strut  depicted  in  Figure  la.  This  strut  is  non- 
•uniform,  the  two  halves  having  different  cross-sectional  areas  and  moments 
thereof. 

The  critical  load  for  a  strut  that  is  unsynnetrical  about  the  center  has 
been  investigated  by  Case, 15  who  has  shown  that  the  value  of  P^  is  given 
by  the  expression 

_2_ _ 1_  (7) 

P  '  ,P  -P 

cr  1  cr  c  cr 

where  iPcr  =  the  crippling  load  of  a  strut  with  two  halves  like  QA,  and 
2Pcr  =  the  crippling  load  of  a  strut  with  two  halves  like  OA  .  It  is  but 
a  matter  of  single  calculation  to  show  that  if  fi  is  the  flexibility  co¬ 
efficient  for  a  centrally  loaded  beam  of  total  length  L  and  of .the  form 
of  OA,  and  if  f2  is  the  value  of  this  parameter  corresponding  to  a  beam 
entirely  like  QA/,  then  the  flexibility  coefficient  for  the  composite  beam 
is  given  by 

2f  =  f.  +  f~  (8) 

x  2 


It  follows  from  equations  (7)  and  (8)  that 


But 


Hence, 


P  f  = 
cr 


,PffT.  ~P  (fn  +  f=) 

1  cr  .  ^  cr  l _ 

,P  +  8P 
1  cr  d  cr 


(9) 


,P  f. 
1  cr  x 


(1C) 


(11) 


It  is  apparent  from  this  that  the  condition  of  complete  uniformity  of 
geometry  is  not  essential.  We  conclude,  too,  that  the  point  of  load  applf 
cation  is  not,  therefore,  the  point  of  symmetry  but  rather  the  point  of 
maximum  compliance^ To  extend  the  argument,  we  examine  next  the  case 
of  a  column  with  continuously  varying  section  geometry.  We  shall  assume 
(for  convenience)  that  the  variation  in  El  is  described  by  the  expression 
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Figure  1 


.  Diagrams  of  Several  of  the  Beam  Configurations 
Studied . 
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11  -  [  1  -  i  <t>2  ]  (12 

and  that  the  end  x  =  0  is  built  in  while  the  end  x  =  L  is  pinned.  It  can 
be  shownl°  that  this  system  has  a  critical  load  that  is  given  by 


17.79 


£  P  £  17.88 
cr 


(13) 


This  strut  has  a  33$  variation  in  the  value  of  El  as  we  proceed  from  one 
end  to  the  other,  a  variation  that  is  greater  than  one  might  expect  in  a 
practical  application. 

Hie  flexibility  coefficient  under  lateral  force  can  be  computed  easily 
from  energy  considerations  (Appendix  i).  There  results 


f  =  0.0109 


and  so 


P  f  =  0.9^5 
cr 


tt^l 


(U) 


(15) 


The  error  between  this  result  and  that  established  in  equation  (6)  amounts 
to  less  than  6$,  which,  when  viewed  from  a  practical  point,  is  excellent. 
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NOHUHIFORM  SUPPORT 

lhe  simple  examples  given  have  provided  substantial  evidence  that  a  practi¬ 
cal  method  for  determination  of  end  fixity  in  columns  exists  and  that  its 
application  would  provide  a  nondestructive  process  for  determining,  in  the 
field,  those  factors  that  are  not  readily  amenable  to  calculation.  But 
there  is  some  weakness  in  the  argument,  since  all  the  results  presented 
apply  to  the  most  mathematically  ideal  boundary  conditions  -  conditions 
that  are  similar  at  both  ends  of  the  member,  a  situation  unlikely  in 
practice.  There  is  lack  of  realism,  too,  in  the  sense  that  we  have  not 
considered  the  variations  in  rotational  and  translational  end  restraint 
which  characterize  virtually  every  frame  structure.  The  analysis  is,  there¬ 
fore,  extended  to  cover  these  situations. 

To  begin,  we  consider  the  case  of  a  uni  form  member  encastre  at  one  end 
and  pinned  at  the  other  (Figure  lb).  This  is  a  classical  problem  and  the 
various  parameters  are  well  established.  The  critical  compressive  load  is 
given  by 


P  =2.05  if  El 


(16) 


The  flexibility  coefficent  at  the  point  of  maximum  compliance  (x  =  O.566) 
is  given  by 


f 


0.009 04 


(17) 


and,  thus,  it  follows  that  the  product 


(ic) 


This  is  not  precisely  of  the  form  of  equation  (fc),  but  the  ■+$  variation 
from  that  value  is,  nevertheless,  small. 

Next,  we  consider  the  situation  shown  in  Figure  1c  of  a  pinned  uniform 
column  but  with  equal  rotational  restraint  at  the  ends.  The  critical 
compressive  load  for  this  case  is  determined  as  a  function  of  spring 
stiffness  from  the  transcendental  equation: 
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:ei 

kL 


tan  u 


(19) 


u 


where 


P 

cr 


El 

12 


(20) 


It  is  obvious  from  the  support  symmetry  that  the  point  of  minimum  lateral 
stiffness  is  the  midpoint  of  the  learn.  The  flexibility  coefficient,  f, 
is  easily  determined,  in  closed  form,  for  this  case  (see  Appendix  II)  and 
is 


1  ~  IfEI  L 


, ,  1EI  x 
i(l  +  “kL  } 


(21) 


The  product  P  f  can  be  written 
cr 


P  f  = 
cr 


2u 


if 


ft  - 


2(1 


2E I,  J 

kl/ 


1^ 


U:) 


It  appears  from  this  result  that  P^f  is  no  longer  a  constant  for  varia¬ 
tions  in  support  but,  instead,  by  reference  to  °quaticn  (•. ).  is  in  error 

by  a  factor 


a 


(21) 


When  the  factor  a  is  plotted  against  spring  stiffness,  k  (Figure  l),  it 
is  clear  that  the  variation  from  unity  is  small  and  that  the  maximum 
error  is  less  than  12$  over  all  values  of  end  fixity. 

These  results  can  be  systematically  extended  to  cases  of  nonsymmetric 
end  restraint.  Suppose  the  previous  example  is  taken  but  with  one  spring 
removed.  The  critical  load  for  this  configuration  is  determined  from  the 
transcendental  equation 


EI  =  I1/_1 
kL  ~  ?u  '  2u 


tan  lu 


(22) 
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Figure 


where  again 


P 

cr 


.  2  El 

4  U  — r 


(25) 


The  flexibility  coefficient  must  now  be  determined  for  given  values  of 
end  restraint.  The  analysis -is  somewhat  more  involved  because  the  point 
of  maximum  value  (minimum  lateral  stiffness)  depends  on  end  restraint.  A 
sunmary  of  the  procedure  used  here  is  given  in  Appendix  III.  The  product 
of  Pcrf  and  the  maximum  flexibility  coefficient  is  plotted  versus  end 
restraint  in  Figure  3.  It  is  readily  apparent  here  that  the  product  is 
approximately  constant  for  unsymmetric  changes  in  end  restraint  and  is 
directly  proportional  to  column  length.  Again,  the  constant  is 

TT^L 

~T5 

and  the  maximum  error  is  less  than  4%  for  all  values  of  restraint.  It  is 
important  from  a  practical  point  to  note  here  that  the  locus  of  points  of 
minimum  stiffness  shown  in  Figure  3  does  not  deviate  from  the  column  mid¬ 
point  by  more  than  10%. 

This  example  can  be  continued  by  taking  the  above  spring  stiffness  as 
infinite  and  then  adding  degrees  of  rotational  restraint  to  the  opposite 
end.  The  critical  load  for  this  case  is  determined  from  the  equation 


?2(u)  =  4*(u)  .  +(u)  ]  (26) 


where 

*l'lL  STaJ  -  31  j  and  *(u)  "  IS  LS  +  ]  (2?) 


and 

P  =  hu2  ^  (28) 

cr  2  v  ' 

L 


Again  the  flexibility  coefficient  and  its  point  of  minimim  value  depend  on 
end  restraint.  When  the  product  Pcrf  is  plotted  against  restraint.  Figure 
4,  it  is  clear  as  before  that  the  result  is  approximately  constant 
throughout  the  complete  range  of  restraint.  The  maxi  mini  deviation  from 
Tl^L A8  is  less  than  4%,  and  the  locus  of  points  minimum  lateral  stiff¬ 
ness  is  within  10%  of  the  column  midpoint. 
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Figure  3*  A  Plot  of  the  Error  Factor,  a 
Against  Rotational  End  Res 
Beam. 
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TTrrERWniATE  LATERAL  SUPPORT  AHD  HIGHER  MODES 


The  consistency  between  the  many  cases  studied  leads  to  the  prospect  that 
a  f  ii'  rnl  but  approximate  relationship  might  exist  between  the  critical 
instability  load  and  the  flexibility  coefficient.  However,  the  fonmxlation 
of  a  general  analytic  procedure  from  which  such  results  can  be  deduced 
appears  somewhat  more  complex  at  present^  therefore,  it  seems  more  effi¬ 
cient  to  study  the  problem  further  by  means  of  specific  exan^les.  Quite 
often,  in  practice,  a  column  will  be  supported  laterally  at  one  or  more 
points,  such  a  configuration  being  adopted  generally  to  in^rove  the  load- 
carrying  ability.  In  some  instances,  the  intermediate  supports  can  be 
considered  as  a  perfectly  pinned  or  clamped  restraint,  in  which  case  the 
column  may  be  subdivided  into  the  intervals  between  supports,  which  in  turn 
may  be  taken  as  columns  with  applied  end  moments.  For  many  other  practical 
cases,  though,  the  intermediate  supports  are  far  from  ideal  and,  for  pur¬ 
poses  of  analysis,  can  be  represented  by  rotational  and  lateral  springs 
located  at  discrete  points.  The  spring  stiffness,  then,  can  be  specified 
to  reflect  closely  the  stiffnesses  of  real  supports.  An  exanqple  of  such 
a  construction  is  depicted  in  Figure  Id,  where  a  pinned  column  is  subjected 
to  a  lateral  spring  restraint  at  its  midpoint.  The  critical  compressive 
axial  load  for  synmetric  deflection  is  determined  when  the  spring  reaction 
tends  to  infinity  and  is  defined  by 

sin  u  !"-  sin  u  +  u(  1  -  l6u2  ^r)  cos  u  1  =  0  (29) 

*-  kLJ  J 


P  =  bv?  ^  (30) 

.  cr  L2 

<*>  (rigid  support),  the  lowest  nontrivial  root  is  u  =  TT, 


where 


Now,  when  k  = 
from  which 


P 

cr  l2 


(31) 


and  when  k  =  0  (no  support),  the  lowest  root  is  u  =  TT/2,  from  which 


„2  El 
P  =  TT  — 5- 
cr  l2 


(32) 


These  are  clearly  the  cases  when  two  waves  and  one  wave  are  produced,  and 
they  establish  limits  on  u.  As  k  is  increased  from  zero,  the  lowest  root 
of  equation  (29)  corresponding  to  symmetric  buckling  is  determined  when 
the  bracketed  term  vanishes.  This  happens  only  so  long  as 


* 
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(33) 


or 


1  -  l6u2  z  0 

kLJ 


kL3  4ff2  El 

16H2  ei  l2 


(3M 


If  this  value  is  larger  than  Pcr  for  a  rigid  support,  then  the  column 
must  buckle  with  two  waves.  The  transition  from  one  to  two  waves  occurs, 
then,  when 


P 


(P  ) 
'  cr 


for  two  waves 


(35) 


or  when 


k  =  16TT 


EI 

t3 


(36) 


For  k  greater  than  this  value,  antisymmetric  buckling  occurs  as  defined  by 
equation  (30).  For  k  less  than  this  value,  symmetric  single-wave  buckling 
occurs  at  a  load  when  (see  equation  (29)) 


2  EX 

-  sin  u  +  u(l  -  l6u  — ^  )  cos  u  =  0 

kLJ 


(37) 


The  solution  can  be  represented  approximately  as 


cr 


a  ♦ 

i6tt  ei  l 


(38) 


kien  the  spring  stiffness  is  zero  or  less  than  that  defined  in  equation  (39), 
the  point  of  minimum  lateral  stiffness  is  clearly  the  midpoint.  The 
flexibility  coefficient  as  shown  in  Appendix  IV  is  then 


f  = 


L3 

48  EI 


[ 


1 


485  EI 


] 


(39) 
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Using  equation  ( _■  ),  there  results 


P  f  = 
or 


m 


where 


a 


j  kL 

1  +  1 6P2  El 


1  + 


kL' 


4£  El 


(41) 


The  factor  a  is  then  a  measure  of  the  deviation  of  Pf  from  the  constant 
value  and  is  plotted  versus  k  in  Figure  5»  It  is  evident  that  the  maxi¬ 
mum  error  is  less  than  7$  and  occurs  at  the  transition  from  the  first  to 
the  second  buckle  mode. 

Mj  2 

Now,  when  r=—  >  16TT  ,  there  is  an  inflection  point  at  l/2,  and  minimum 

lateral  stiffness  clearly  does  not  occur  here.  It  does  occur,  however,  at 
two  points  located  midway  between  the  ends  and  center  support.  This 
requires  some  modification  of  our  previous  results.  Consider  the  simulta¬ 
neous  application  of  a  point  load,  Q,  at  one  of  the  points  and  an  equal 
and  oppositely  directed  load  at  the  other.  The  sum  of  the  resulting  de¬ 
flections  is 


-L 


#EI 


(42) 


and,  if  we  define  a  new  flexibility  coefficient 


then 


192  El 


1*  tT~  ei 
.2 


(43) 


(44) 


192  EI 


which  is  just  the  same  result  as  that  for  the  earlier  examples 


This  approach  can  also  be  extended  to  other  situations  in  which  higher 
modes  are  critical.  The  first  example  taken  is  that  of  a  column  that  is 
encastre  at  both  ends  and  pinned  at  a  center  support.  In  this  case,  the 
buckle  mode  is  again  antisymnetric  and,  using  equations  (l6)  and  (17), 
the  points  of  minimum  stiffness  occur  at  X  =  0.29  L  and  0.71  L.  If 
oppositely  directed  loads,  Q,  are  applied  at  these  points,  then  it  follows 
in  a  similar  manner  to  the  previous  case: 


6  =  6  + 


2  (.00984) 


(45) 


I 


and 


6 

Q 


0.00246 


(46) 


The  critical  compressive  load  is 


P  =  2.05TJ2 
cr 


El 

(|)2 


=  8.2T12  ^ 
L 


(47) 


and  so 


P  f0 
cr  2 


97  A 


(48) 


giving  an  error  of  3$* 

If  we  consider  buckling  in  the  third  mode  for  a  pinned-ended  column,  then 
one  can  show  that  three  alternating  forces  are  required  and  that  if 


th^n 


f3 


6  +6+6 
12  3 


(49) 


V.. 


p  f 


cr 


T^L 

~TE 


(50) 


with  zero  resulting  error. 
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It  is  possible  to  treat  higher  mode  buckling  cases  with  the  various 
boundary  conditions  discussed  earlier  by  subdividing  the  column  at  inflec¬ 
tion  points  into  a  series  of  column  elements.  As  before,  the  flexibility 
coefficient  must  be  defined  by  applying  forces,  to  each  element  at  points 
of  minimum  stiffness  and  summing  the  magnitudes  of  the  deflections  pro¬ 
duced;  then 


f 


n 


(5D 


and 


P  f 


cr 


TT^L 


(52) 


The  Pf  product  is  exact  for  pihned-ended  conditions.  Since  for  higher 
mode  buckling  with  other  end  conditions,  the  interior  portion  of  the 
column  deforms  into  sinusoidal  waves,  it  can  be  concluded  that  the  error 
for  these  cases  diminishes  with  increasing  buckle  mode. 

These  results  suggest  a  generalized  procedure  for  situations  when  the 
buckling  load  corresponds-to  higher  modes.  The  maximum  generalized  flexi¬ 
bility  coefficient,  rn  =  X  6i/Q,  mist  be  determined  forkall  contemplated 
modes  by  first  applying  a  single  concentrated  lateral  load  and  then  adding 
additional  loads  in  an  alternating  fashion,  in  each  case  adjusting  the 
positions  of  the  loads  so  that  the  deflection  sum  is  a  maximum.  This  de¬ 
fines  the  "maximum  flexibility"  configuration  for  that  mode,  and  the  mode 
with  the  lowest  value  will  then  be  preferred,  so  that 


(53) 


P  f 
cr  n 


~  Tfx 
=  TH 
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f  LATERAL  RESTRAIHT  AT  OHE  EHD 

We  have  illustrated  by  the  several  examples  that  when  the  end  fixity  condi¬ 
tions  as  regards  rotational  restraint  are  varied,  the  product  of  the 
critical  load  and  the  influence  coefficient  at  the  point  of  wnrtmim  compli¬ 
ance  is  very  nearly  constant  for  both  uniform  and  nonuniform  columns.  We 
now  turn  to  Investigate  the  influence  of  variation  in  lateral  restraint  at 
one  end.  To  study  this,  ve  consider  the  column  that  is  built  in  at  its 
root  but  is  free  at  the  tip.  The  critical  joad  for  such  a  column  when  El 
is  constant  along  the  length  is  given  by 


P 

cr 


(5*0 


The  point  of  maxlmim  compliance  is  clearly  the  tip  of  the  column,  and  the 
flexibility  coefficient  is  from  the  well-established  deflection  formula 
given  by 


f  = 


(55) 


Thus,  the  product 


P  f  = 


cr 


L3  _  TfT, 

3  El  ~  12  El 


(56) 


There  is  clearly  a  large  discrepancy  here.  It  is  interesting,  therefore, 
to  determine  at  what  point  along  the  column  the  flexibility  coefficient 
must  be  determined  for  this  product  to  have  the  same  value  as  for  the 
other  cases  considered.  It  follows  that,  since  the  coefficient  varies  with 
the  cube  of  the  distance  of  the  load  point  from  the  root  of  the  column,  the 
corresponding  point  is  defined  by  the  length  L^,  where 

Lx  =  ~  .625L  (57) 


It  is,  perhaps,  worthwhile  to  note  that  this  point  is  not  far  removed  from 
that  which  is  determined  for  the  case  when  the  tip  is  pinned.  In  that 
case,  the  point  of  maximum  compliance  was  located  at  .586  L. 

To  what  degree  this  near  correspondence  would  be  influenced  by  a  reason¬ 
able  nominiformity  of  section  is  now  studied.  The  case  considered  here  is 
that  of  a  cantilevered  column  in  which  the  cross-sectional  characteristic 
is  refined  by 


1 


A. 


20 


(58) 


El  =  El 


1  5  L2  ^ 


where  EI0  is  the  value  of  the  section  parameter  at  the  clawed  end  and  x 
is  measured  from  that  point.  The  critical  load  for  this  structure  is  given 


by 
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P  = 
cr 


12  El 


5L 


(59) 


The  deflection  due  to  a  unit  load  (the  flexibility  coefficient)  positioned 
at  a  point  XL  from  the  fixed  end  is  given  by 


f  = 


XL  2 

r  fo-d — r  <u 

•  no  [ 1  -  ?  <!>  ] 


(6o) 


The  point  is  to  determine  the  value  of  X  which  makes  the  product  Pcrf  have 
the  value  previously  established.  For  this  condition  to  be  met,  the  value 
of  f  must  be  given  by 


f  -  5TT2  .  l£_ 
f  (12)  (to)  EIq 


(61) 


A  solution  of  equations  (59)  and  (60)  yields  the  value 


X  =  0.634  (62) 

or 

L1  =  0.634L  (63) 

which  is  remarkably  close  to  the  value  0.625  obtained  for  the  uniform 
cantilever,  equation  (56). 
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I*  appears,  then,  for  cantilever- type  configurations  of  uniform  or  non- 
uniform  geometry,  that  the  product  Fcrf  can  be  taken  as  constant  with 
reasonable  error,  provided  the  flexibility  coefficient  is  measured  within  _ 
13$  of  the  column  midpoint.  When  the  maxi  mini  value  of  the  coefficient 
falls  outside  this  limit,  then,  the  value  at  the  limit  must  be  used  Instead. 
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CONCLUSIONS 


The  research  outlined  in  this  paper  has  established  a  practical  yet  simple 
method  of  determining  the  end  fixity  of  realistic  columns.  It  has  been 
shown,  both  analytically  and  experimentally,  that  the  end  fixity  coeffi¬ 
cient  for  a  column  can  be  obtained  within  usual  engineering  tolerance  by 
dividing  the  parameter  TT^l/^B  by  the  product  of  the  Euler  load  and  the 
maxi  mini  value  of  the  flexibility  coefficient  for  a  point  lateral  load 
applied  between  3/8  and  5/8  of  the  span  of  the  column.  This  rule  is  appli- 
cable  for  all  conditions  of  restraint,  at  either  end  of  the  column,  so 
long  as  the  applied  force  is  constrained  to  remain  vertical. 
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APPENDIX  I 

FLEXIBILITY  COEFFICIENT  FOR  A  TAPEREL  BEAM 


Consider  the  clamped -pinned  beam  shown  in  Figure  lb. 
El  is  variable  and  is  given  by 

4 

El  =  El  (l-a2^2) 
o  5 


The  section  stiffness 


(64) 


where  %  =  x/L  and  where  a  is  the  taper  ratio.  It  is  desired  to  find  the 
position  X  for  which  the  flexibility  coefficient  associated  with  loads 
and  displacements  at  x  =  XL  is  maximized.  The  p~  oedure  employed  here  is 
the  well-known  d ninny  load  method. 

The  first  requisite  in  the  analysis  is  to  obtain  the  actual  moment  distri¬ 
bution,  for  which  it  is  necessary  to  determine  the  reaction  force,  R,  at 
the  end  of  the  beam.  Toward  this  objective,  the  internal  complementary 
work  is  given  by 


* 

U  = 


L 

r 

j 

o 


(65) 


f 

/ 

Since  the  reactions  do  no  work,  the  complementary  potential  is  given  by 
an  expression  identical  to  equation  (65).  In  terms  of  the  unknown  redun¬ 
dant  R,  the  moment  distribution  is  given  by 


rRL  (1-S)  -  QL  (X-5) 
^  =  iRL  (1-5) 


0  *  5  £  X 
X  £  §  £  1 


(66) 


Substituting  equation  (66)  and  equation  (64)  into  equation  (65)  gives, 
after  rearrangement. 


* 

U  ^ 


o  i-a  r 


(67) 
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Now,  the  force  at  the  end  is  that  R,  among  all  R's  in  equilibrium  with  the 
later  load  Q,  which  permits  compatible  deformations  of  the  structure; 
equivalently. 


Carrying  out  the  indicated  differentiation  of  equation  (67)  and  solving 
for  R,  we  have 


(R/Q)  = 


>af  gs 

l-a2§ 


The  deflection,  6,  of  the  beam  is  given  by  the  dummy  load  method  as 


r  W* 


where  Mq  is  the  moment  distribution  due  to  the  dumny  load.  In  this  case, 
the  dummy  load  is  applied  at  the  same  place  as  the  actual  load,  Q,  so 


Mo  (*)  =  ^  M1  W 


Substitution  into  equation  (70)  gives,  after  rearrangement. 


(6/Q)  = 
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/ 

Thus 


Now,  equation  (69)  can  be  integrated  in  closed  form  to  yield 
R,(R/q)=  a~\  In  (~  +  2a  (1  +  \)  In  VI  -  a2\2 

+  2  In  .  \mL -  2  a\  (7*0 

VI  -  a“\~ 


where 

Rx  =  a2  In  f  4a  Id  VT -  a2  +  2  In  ^  *  a  ij  -  2a  (75) 

Furthermo substitution  of  equations  (7*4)  and  (75)  into  equation 
(73)  gives,  after  considerable  simplification. 


■  3  L7q)_ 

(LJ/EI  ) 


2  2 

a  \ 


In 


1  +  a\ 
1  -  a\ 


+■  4a\  In 


~Z7Z 

a  \ 


+  2  In 


1  +  a\ 

/i  -  .V 


2a K  ( 76) 


Since  it  was  desired  to  maximize  the  influence  coefficient,  it  was 
necessary  to  differentiate  equations  (74)  and  (76).  Carrying  out  the 
differentiation  yields,  after  simplification. 


5 

n  ’  rf  v 


(7q)  =  a-  m  ( 


1  +  a\ 
1  -  a\ 


)  +  2a  In 


r.  272 

V  1-a  \ 


(77) 
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or 


a 


2  A 

d\ 


- 3  l>-  (E/«  ]  *  <i^£  > 


(L  /EIq) 


2  [l  -  (R/Q)  ]  In  \l  1  -  aV 


(78) 


Hence,  the  problem  of  maximizing  the  influence  coefficient  is  expressed  by 


f  W-*2K  t^T  J  ■  0  <'9) 

IL  /EI  ; 


Again,  it  is  apparent  that  the  root  cannot  be  found  in  closed  form;  a 
numerical  scheme  is  necessary. 

Newton's  method  was  used  here,  with  the  iteration  being  accomplished  via 
a  digital  computer.  The  maximum  flexibility  coefficient  and  the  associated 
position  for  a  taper  ratio,  a,  of  l/3  were  found  to  be: 


\  =  0.599 


(80) 


f  =  1.09  x  10 


(81) 


Taper 

Maximizing 

Pinned-End 

Flexibility 

Ratio 

Lambda 

Reaction 

Coefficient 

T 

.589456 

.416455 

1.01000 

.2 

.593348 

.418838 

1.04138 

•  3 

.597658 

.421547 

1.07608 

.4 

.602422 

.4246 

1.11463 

•5 

.607749 

.428089 

1.15850 

.6 

.613806 

.432163 

1.20853 

.7 

.620837 

.437029 

1.26699 

.8 

.629294 

.443103 

1.33743 

.9 

.640046 

.451192 

1.42667 
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APPENDIX  II 


FLEXIBILITY  COEFFICIENT  FOR  A  BEAM  WITH  EQUAL 
_ ROTATIONAL  END  RESTRAINT 


The  problem  considered  is  to  determine  the  maximum  flexibility  coefficient 
of  a  pinned  beam  with  rotational  restraint  as  shown  in  Figure  lc.  A  com¬ 
plementary  energy  approach  similar  to  that  of  Appendix  I  was  used.  The 
problem  is  symmetric  about  the  beam  midpoint,  which  is  clearly  the  poino 
of  maxi  naan  flexibility  coefficient.  The  moment  distribution  can  be  written 


M  =  M  + 


& 

2 


0  £  x 


£  L/2 


(82) 


so  the  internal  complementary  energy  is 


* 

U  = 


L 

Bei 


+  M  QL  + 

o 


(83) 


The  work  done  by  the  spring  reaction  forces  is  -  >^0,  where  0  is  the  end 
rotation.  The  total  complementary  potential  is  then  the  sum 


TT*  =  U*  -  M  0 
o 


(84) 


The  proper  value  of  Mq  acting  in  the  structure  is  determined  by  the  prin¬ 
ciple  of  minimum  total  complementary  potential: 


0  = 


(85) 


Using  equation  (84)  in  equation  (85)  and  solving  for  Mq  yields 


M 

o 


(86) 


The  midpoint  deflection  can  now  easily  be  found  either  by  the  dunny  load 
method  or  equivalently  by  Castigliano'r  First  Theorem: 


V, 

’t 


Using  the  latter  approach  requires  determination  of  the  strain  energy  U: 

i 

2 


u  -  rs  (""o2  +  Mo«L  ♦  +  % 


2k 


(88) 


Now,  using  equations  (86)  and  (88)  in  equation  (87)  gives,  after  some 

manipulation. 


6 "  ufe  ( 1 '  Kl  i  J 

J  v  kL' 


(89) 


at  the  beam  midpoint,  from  which  the  flexibility  coefficient  is  determined 
as 


f-i-  L' 


q  Wei 


[  1  "  .  EI\  ] 


^  +  E> 


(90) 


( 


APPENDIX  III 

FLEXIBILITY  COEFFICIENT  FOR  A  PINNED  BEAM  WITH  ROTATIONAL 

RESTRAINT  AT  ONE  END 


The  maximum  flexibility  coefficient  for  the  pin-ended  beam  with  rotational 
restraint  at  one  end  (Figure  lc  with  the  leftmost  spiring  removed)  can  be 
determined  with  good  accuracy  by  a  Rayleigh-Ritz  procedure.  This  method 
consists  of  choosing  an  approximate  ninction  for  the  deflected  shape  of 
the  beam  and  then  adjusting  the  unspecified  parameters  in  the  function  by 
minmizing  the  total  potential  energy  in  the  system. 

For  the  present  case,  ■'.he  deflection  function 


N 

E 

m=l 


»  .  mTlx 

A  sin  -=— 
m  L 


(9D 


consisting  of  the  first  N  sin  terms  was  chosen.  The  total  potential  energy 
can  t^en  be  written 


TT  =  M 


i1, 


A 

a/ 


)  dx 


(92) 


The  proper  values  of  the  constants,  are  determined  by  requiring  that 
iT  assume  a  minimum  value;  i.e.. 


for  i 


1, 2  • * . «N 


(93) 


Substituting  equation  (9l)  in  equation  (92)  and  carrying  out  the  indicated 

operations  yields  an  algebraic  equation.  Using  this  in  equation  (93),  then, 

gives  a  set  of  simultaneous  linear  equations  determining  the  constants, 

A  ,  as  follows: 
n 


(x) 


A 

m 


nh  cos 


sin  (\mIT)  for  m  =  1,2  ....N 


(94) 
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These  equations  can  be  solved  in  a  straightforward  manner  for  the  constants. 
Am,  given  a  valne,  for  the  load  position.  For  the  present  case,  a 
direct  numerical  solution  was  c 01153 le ted  and  the  resulting  coefficients,  V 
were  used  to  find  the  flexibility  coefficient  at  the  point  \L: 


6  _ 

Q  Q 


(95) 


Values  of  \  were  adjusted  in  an  iteration  scheme  until  the  flexibility 
coefficient  assumed  a  maximum  value.  The  results  are  presented  in  Figure  3* 
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APPENDIX  IV 

FLEXXBILI1Y  COEFFICIENT  FOR  A  PINNED  BEAM 
WITH  INTERMEDIATE  LATERAL  SUPPORT 


The  maximum  flexibility  coefficient  for  the  pinned  beam  shown  ir.  Figure  Id 
with  a  lateral  spring  support  at  its  midpoint  can  be  determined  In  a  direct 
manner  by  energy  methods.  For  values  of  spring  stiffness  less  than  that 
given  in  equation  (3 6),  i.  e.. 


k  < 


161^ 


(96) 


the  point  of  maximum  flexibility  coefficient  (minimum  stiffness)  is  the 
midpoint.  If  a  load,  Q,  is  applied  here,  the  moment  distribution  is 


M  =  Q  ~  R  x  0  s  x  s  L/2 

=  x  -  (Q  -  R)(x  -  L/2)  §  *  x  *  L  (97) 


where  R  is  the  upward  reaction  due  to  the  spring.  Essentially,  the  same 
energy  procedure  outlined  In  Appendixes  I  and  II  is  used  here.  The  total 
complementary  potential  can  be  written 

l4 

TT*  =  eJ  j*  x)2  dx  -  R(-6)  ‘  (98) 

b 


where  6  is  the  midpoint  deflection  *n  the  downward  direction  of  Q.  As 
before,  the  unknown  reaction,  R,  has  that  value  which  minimizes  IT ;  equiva¬ 
lently, 


0  = 


(99) 


Using  equation  (96)  in  equation  (99)  and  simplifying  gives  the  result 


R  = 


1  + 


k _ 

kd  El 

kL3 


3^ 


(100) 


The  deflection  of  the  spring  is  just  the  deflection  at  the  midpoint  due  to 
Q.  >  Thus, 


■I  t-^3-] 


1  4. 

1  bb  El 


(101) 


and  the  flexibility  coefficient  is  given  by 

"  l3  r  1  ' 

kL3  - 

1  TiSlrf 


(102) 
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